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ABSTRACT 

The scale radius of dark matter halos is a critical parameter for specifying the density distribution of dark 
matter, and is therefore a fundamental parameter for modeling galaxies. We develop here a novel, observation- 
ally motivated probe to quantitatively infer its value. We demonstrate that disturbances in the extended atomic 
hydrogen gas disks of galaxies can be used to infer the scale radius of dark matter halos. Our primary metric 
is the phase of the m = 1 mode of the disturbance in the outskirts of the gas disk, which we take to be produced 
by a tidal interaction. We apply the method to the Whirlpool Galaxy, which has an optically visible satellite. 
We explore potential degeneracies due to orbital incUnation and initial conditions and find our results to be 
relatively insensitive to these considerations. Our method of tracing the dark potential well through observed 
disturbances in outer gas disks is complementary to gravitational lensing. 

Subject headings: Galaxies: evolution - Galaxies: dynamics, Dark matter - indirect methods 



1. INTRODUCTION 

The cold dark matter (CDM) paradigm of structure for- 
mation is successful at recovering the basic skeletal struc- 
ture of the universe - the large-scale distribution of galaxies 
(Springel et al. 2006). However, the agreement between the- 
ory and observation is less secure when this model is applied 
to galactic (and sub-galactic) scales. While the existence of 
dark matter halos in galaxies was observationally inferred in 
the 1970's (Rubin et al. 1977; Bosma 1978), there are few 
observational probes that can be used to infer the details of 
the distribution of dark matter. Simulations have not yet been 
able to resolve the mitigating effects of gas cooling (Gnedin 
et al. 2004) versus outflows (Govemato et al. 2010), which 
can act to increase (decrease) the distribution of dark matter 
in galaxies. In the coming decades, the CDM paradigm will 
be increasingly vetted against detailed observations of the in- 
nards of galaxies that will delve even more deeply into galaxy 
evolution. 

Navarro, Frenk & White (1996; 1997; NFW) found that 
dark matter halos formed through dissipationless hierarchi- 
cal clustering have an universal density profile. Due to the 
absence of scale in this heirarchical structure formation sce- 
nario, dark matter haloes are scaled versions of each other, 
with the total mass being the scaling parameter. The NFW 
density profile is given by: 



(r/R,){l+r/R,)^ 

where Pc is the critical density for closure of the universe, and 
Sc is a characteristic overdensity for the halo. The scale radius, 
Rs = Rioo/c is a characteristic radius, and c is the concentra- 
tion parameter that relates this characteristic radius and Rzm, 
which is the radius at which the mean enclosed dark matter 
density is 200 times the critical density. The characteristic 
overdensity of the halo, Sc, is simply a function of the con- 
centration parameter c (NFW). The NFW profile differs from 
the Hernquist (1990) profile in its asympotic behavior, scal- 
ing as r"^ rather than r""* for r/r^ ^ 1. The scale radius of 
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the dark matter halo along with the concentration parameter 
specify the density distribution for a given halo mass in the 
NFW formalism. Currently, these parameters are prescribed 
from large cosmological N-body simulations (Bullock et al. 
2001; Maccio et al. 2008). It is not presently clear how we 
can determine these values for individual galaxies. What we 
seek to do here is to develop an independent, observationally 
motivated probe of the scale radius of dark matter halos, and 
apply it to the well-known Whirlpool Galaxy (M51) to in- 
fer its scale radius. The density switches its behavior at the 
scale radius, from r"' for r ^ R^ to for r ^ R,, and passes 
through the transitional region for r ^ R,. Thus, the scale 
radius of dark matter halos is a critical parameter for speci- 
fying the density distribution of dark matter in galaxies, and 
therefore a fundamental parameter for modeling galaxies. 

The extended atomic hydrogen (HI) disks of galaxies pro- 
vide an unique probe of galaxy evolution. They are ideal 
tracers of tidal interactions with satellites and the galactic 
gravitational potential well. We recently developed a novel 
method whereby one can infer the mass, and relative position 
(in radius and azimuth) of satelUtes from analysis of observed 
disturbances in outer gas disks, without requiring knowl- 
edge of their optical light (Chakrabarti & Blitz 2009, hence- 
forth CB09; Chakrabarti & Blitz 2011, henceforth CBll; 
Chakrabarti, Bigiel, Chang & Blitz 2011, henceforth CBCB; 
Chang & Chakrabarti 2011, henceforth CCll). We apphed 
this method to M5 1 and inferred that its satellite has a mass 
one-third that of the primary galaxy, with a pericentric ap- 
proach distance of 15 kpc. We found these estimates to be 
corroborated by observations (Smith et al. 1990) and recent 
simulation studies (Salo & Laurikainen 2000; Dobbs et al. 
2010). Moreover, at the time when our simulations achieve 
the best-fit to the HI data, the azimuth of the satelhte in the 
simulations agrees closely with its observed location. CBCB 
note that the derivation of these numbers is uncertain at the 
factor of two level due to variations in the initial conditions 
of the simulated M51 galaxy, orbital inclination and orbital 
velocity of the satelhte. We call this method "Tidal Analysis". 

We frame the problem here in the following way - for 
galaxies like M51, for which we have demonstrated that our 
Tidal Analysis method works, can we determine the scale ra- 
dius of the dark matter halo of M51? We take the satelhte 
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mass and pericentric distance from our earlier studies as an in- 
put. It is worth noting that such a procedure has viable and im- 
mediate applicability. The THINGS survey produced HI maps 
of local spiral galaxies (Walter et al. 2008), many of which 
display disturbances in the outskirts (Bigiel et al. 2010). At- 
tempting to constrain the dark matter halo in this way is sim- 
ilar in spirit to earlier work on analysis of stellar tidal tails by 
Mihos, Dubinski & Hemquist (1998). If the satellite mass and 
pericentric distance of the satellites that produce these distur- 
bances can be characterized, as for M51 (CBCB), then as we 
show below, the phase of the m = 1 mode directly yields the 
scale radius of the dark matter halo. This method of probing 
the dark matter mass distribution is independent of the stel- 
lar light and is complementary to strong gravitational lensing, 
which primarily probes regions interior to the Einstein radius, 
^10 kpc in spirals (Wright & Brainerd 2000; Treu & Koop- 
mans 2002). The weak lensing signal has also been exploited 
by stacking the surface density contrast to produce mean den- 
sity profiles for clusters (Sheldon et al. 2008) and total galaxy 
masses (Mandelbaum et al. 2006). This paper is organized 
as follows: in §2, we demonstrate that the phase of the m = 1 
mode can be used to quantitatively infer the scale radius of 
dark matter halos and apply the method to M5 1 . We discuss 
caveats and future work in §3, and conclude in §4. 

2. RESULTS 

The simulations we discuss here have the same setup as in 
CBCB. Since that paper discusses the simulation setup in de- 
tail, we only briefly discuss it here. We carry out SPH simu- 
lations using the GADGET code (Springel 2005), of M51 in- 
teracting with its companion. CBCB carried out a simulation 
parameter survey and compared the resultant Fourier ampli- 
tudes of the low order modes of the gas surface density with 
the observed HI data of M5 1 . They found that placing sim- 
ulations on a variance vs variance plot (where the variance is 
with respect to the low order modes of the simulations and the 
data) made the best-fit simulations visually apparent. Specif- 
ically, CBCB found that the best-fit to the HI data occurred 
for a 1:3 mass ratio satellite with a pericentric distance of 
15 kpc, parameters that are corroborated observationally and 
are in agreement with other simulation studies (Smith et al. 
1990; Dobbs et al. 2010; Salo & Laurikainen 2000). CBCB 
also found that the azimuthal location of the companion of 
the best-fit simulation agrees very closely with the observed 
location of M5rs companion, at the time when the Fourier 
amplitudes most closely match the data. Thus, CBCB con- 
cluded that analysis of observed disturbances in the extended 
HI disks of galaxies can be used to infer the mass and current 
distance (in radius and azimuth) of galactic satellites. Ear- 
lier work in this series of papers presented the basic reasoning 
as to why the mass-pericentric approach degeneracy in the 
tidal force can be broken when the time integrated response 
of the primary galaxy is considered (CB09), and described 
the method to find the azimuth of galactic satellites from the 
phase of the modes (CBl 1). 

Our main goal in this paper is to determine whether the 
scale radius of the dark matter halo in the primary galaxy can 
be inferred from analysis of observed disturbances in the ex- 
tended HI disk of M51. Since we have earlier characterized 
the mass and pericentric approach distance of M5 1 's compan- 
ion, we take these quantities as inputs in our study here. Here, 
we primarily vary the density profile (and hence the potential 
depth) of the dark matter halo of M5 1 , which we take to fol- 
low an NEW profile. We investigate whether varying the scale 



radius of the dark matter halo of M5 1 will be reflected in the 
disturbances of the extended HI disk cleanly enough to allow 
us to infer its value. 

We begin by studying the gas density response as M5 1 in- 
teracts with a 1:3 mass ratio satellite with a pericentric dis- 
tance of 15 kpc (parameters we derived in CBCB), while we 
vary the scale radius of the dark matter halo. We vary the scale 
radius from low to high values {R, = 1 1 - 32 kpc) to investi- 
gate its effect on the resultant disturbances in the HI disk. The 
scale radius is related to the concentration parameter and the 
outer radius of the dark matter halo, i.e., = Roao/c. There- 
fore, we can either hold /?200 constant and vary the concentra- 
tion parameter, or hold the concentration parameter constant 
and vary /J200- The mass of the dark matter halo scales as R\qq 

(M200 = 200pe j7r7?2oo)' while the concentration parameter is 
related to the angular momentum of the halo as motivated 
by the Mo, Mao & White (1998) formalism (Springel et al. 
2005), and affects the size of the baryonic disk. Our earlier 
models (CBCB) were based on the /?20{) =160 h"' kpc case, 
which gives a galaxy mass consistent with observational esti- 
mates (Leroy et al. 2008). We first set Rom = 160 h~' kpc and 
vary the concentration parameter, which is equivalent to vary- 
ing the scale radius of the dark matter halo. Below we show 
that varying the concentration parameter and R200 such that 
the scale radius is constant gives nearly identical results for 
the phase of the m = 1 mode in the outskirts, which demon- 
strates that the critical parameter that governs the formation 
of these disturbances (once we have an observational handle 
on the mass of the galaxy from the rotation curve) is the scale 
radius. Thus, the two cases we focus on here are: 1) holding 
/?2()() constant and varying c which corresponds to varying R/, 
this means we hold the mass of the dark halo constant while 
we vary the scale radius, where the density switches from r"' 
for r ^ Rs to for r ^ Rs, and 2) holding Rs constant and 
varying /?200 (which will then also vary c) and therefore the 
mass of the dark halo. 

Figure [T] shows the resultant gas density response of M5 1 
when we vary the scale radius (concentration) of the dark mat- 
ter halo of M51 from a small (large) value (R^ = ll,c = 14), 
to the fiducial value used earlier by CBCB {R, = 17, c = 9.4), 
to a low (high) value (Rs = 32, c = 5). Varying the scale radius 
varies the potential depth of M5 1 and will therefore dramati- 
cally impact the formation of tidal features, as we see clearly 
from Figure [1] Steeper potential wells, which are produced 
by steeper density profiles, are more effective at holding on to 
the gas (Mihos, Dubinski & Hernquist 1998). Therefore, the 
Rs= 1 1 , c = 14 case yields more tightly wound structures rel- 
ative to the fiducial value of the scale radius (R, = 17, c = 9.4), 
as well as of course the largest scale radius we consider here 
(Rs = 32, c = 5), where the density profile follows the shallow 
slope of r"' nearly all the way out to where we can probe the 
extended HI disk of M5 1 . 

Figure |2|a) shows how the phase of the m = I mode varies 
as we vary the scale radius of the dark matter halo in the three 
cases shown in Figure [T] For projected gas surface density 
denoted S(r, 0), we calculate the phase of individual modes 
"m" by taking the EFT, where below we have set m= 1: 



(j)(r, m) = arctan 



[-ImagFFT ^(r^cj))] 
[-Re FFT E(r, cj))] 



(2) 



The phase of the modes contains information on the shape 
of the spiral planform, i.e., tightly wrapped spirals will have 
a sharp gradient in the phase, while open spirals will have 
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a flatter profile (Shu 1983; Chakrabarti & Blitz 2011). For 
the fiducial value of = 17(c = 9.4) used by CBCB (shown 
in the solid green line), the phase of the m= I mode in the 
simulation is fairly flat, as is that of the HI data of M5 1 , shown 
in the dash-dotted green line. However, for the other cases, 
i.e., 7?j = 11 and = 32, the phase is either steeply falling or 
rising. For the Rs = 32 kpc case, the density profile follows 
a shallow r"' slope nearly all the way to the edge of the HI 
disk of M51, and as such this is the loosest spiral (nearly a 
straight line) produced. The other extreme is the /?, = 1 1 kpc 
case where the density follows a profile in the radial range 
shown in Figure|2|a). Since the density profile is much steeper 
for the latter case, it is harder to "unwrap" the spiral, leading 
to a tighter spiral planform and positive gradient in the phase 



Fig . 1 . — A series of controlled experiments showing how M5 1 's gas den- 
sity response varies as the scale radius of the dark matter halo is varied. We 
hold R200 constant here as we vary the concentration parameter (scale radius) 
from large (low) to high (small) values. The mass of the satellite and peri- 
centric approach are taken from CBCB. From top to bottom the cases are: (a) 
= 1 1, c= 14, (b) Rs = 17, c = 9.4, the standard case, and (c) Rj = 32, c = 5. 
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Fig. 2. — (a) Phase of the m = 1 mode as the scale radius of the dark matter 
halo is varied, shown at the same times as the gas density images. The three 
cases correspond to probing regions inside the scale radius (the low concen- 
tration case), close to the scale radius (the standard C=9.4 case), and outside 
the scale radius (a high concentration case), (b) Phase of the m = 1 mode as 
a function of r/R, for the S, = 17 case. Marked are the three regimes of in- 
terest. d/dr{4>i ) < Ofor r/S., < 1, transitions at r/Rs ~ 1, and d/dr(4>i ) > 
for r/Rs > 1 until the edge of the gas disk. R200 is held fixed here as the 
concentration is varied. 



of the m = 1 mode. The radial variation of the phase, therefore, 
gives us a handle on the scale radius. 

There are three radial regions of interest in the behavior of 
the phase, and we mark them explicitly in Figure [2|b) for the 
Rs = 17 kpc case. Here, we plot the phase as a function of 
the dimensionless variable r/Rs- We see that d/dr{(f)i) < 
when r/Rs < 1, there is a transition region for r/Rs ^ 1, and 
d/dr{(f)i) > when r/Rs > 1, until the edge of the gas disk 
is reached. In §2.1, we show that the (slight) uncertainty in 
the phase that arises due to variations in initial conditions and 
orbits does not preclude us from differentiating between the 
cases shown in Figure |2] If one can observe the transitional 
region, where the gradient in the phase switches from declin- 
ing to increasing, then one can infer the value of the scale 
radius. The shape of the phase of the Rs= 17 kpc case is quite 
flat, as is that of the HI data for M51. This comparison there- 
fore indicates that the scale radius of the dark matter halo in 
the Whirlpool Galaxy is ^ 17 kpc. 

Our inference of a scale radius of 17 kpc for M51 is com- 
parable to values expected from dissipationless cosmologi- 
cal simulations (Bullock et al. 2001). These simulations 
predict a 68 % spread in concentration values from 9.2-21.3 
for virial mass halos of ^ lO'^M© having virial radii /?vir ^ 
200 h"' kpc, which would yield Rs ^ 13-31 kpc. More re- 
cent work (Maccio et al. 2008) finds similar concentration- 
mass relations. The NFW profile emerges in simulations that 
do not include baryonic physics. One may well expect that 
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baryonic processes, such as gas cooling or supernova feed- 
back (Gnedin et al. 2004; Governato et al. 2010) would alter 
the density profile of the dark matter halo. It is reassuring 
therefore that derivation of the scale radius for a specific lo- 
cal spiral with this new method (that relies on an analysis of 
disturbances in the extended gas disk) yields a value that is 
consistent with the expected range from dissipationless, cos- 
mological simulations. 



2.1. Dependence on Initial Conditions 
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Fig. 3 . — The variation of the phase as the concentration is varied (by vary- 
ing /?20o) while holding the scale radius of the dark matter halo constant. 
The black symbols mark the C = 9.4, Rj = 17 case, the orange marks the 
C = 10.5, R, = 17, and the red the C = 11.7,Rs = 17. As is clear, the be- 
havior of all three is nearly the same for r/S, > 1, and all three have similar 
slopes in the three regimes of interest noted in Figure |2] For r/R, < 1, the 
different concentration parameters yield different disk sizes, which is why the 
amplitude of the phase is different in the inner regions. 



It is important to emphasize that for r/R^ > 1, our results 
depend essentially on the scale radius of the dark matter halo, 
and not on the concentration parameter. We demonstrate this 
in Figure[3] where we plot the phase of the m=l mode of three 
cases where we vary the concentration parameter and R200, 
but hold the scale radius fixed. As is clear, for r/Rs > 1, the 
phase depends primarily on the scale radius (until we reach 
the edge of the gas disk; the size of the gas disk does depend 
on the concentration parameter, which is why there are small 
differences at the very edge). All three cases, for fixed R,, 
have the same type of slope in the three regimes of interest, 
i.e., for r/Rs < 1, d/dr{(f)\) < 0, there is a transition region 
for r/Rs ^ 1, andd/dr((t)i) > when r/R, > 1. For r/R, < 1, 
the concentration parameter affects the size of the disk, and 
hence dramatically affects the amplitude of the phase in the 
inner regions. 

Ideally, one would observe the transitional region where 
the phase of the m = I mode switches from d/dr{(f)\) < for 
r/Rs < 1 to d/dr{(t>\) > for r//?j > 1 to infer the value of the 
scale radius for a specific tidally interacting galaxy. If R, is so 
large that it is comparable to the extent of the HI disk (Rhi), 
then one can at best put a lower limit on the scale radius, i.e., 
if we observe that the phase of the m = 1 mode continues to 
have a negative gradient out to Rhi, then R, is greater than 
Rhi- hi disks are generally more extended than the optical 
radius (Wong & Blitz 2002), so the simulation modeling of 
the baryonic regions can in general be informed by optical 
or near-infrared photometry to yield stellar masses and disk 
sizes. 




Fig. 4. — Phase of the m = 1 mode as the scale radius of the dark matter 
halo is varied from (a) R, = 1 1, c=14, (b) Rj = 17, c = 9.4, the standard case, 
and (c) R, = 32, c= 5, shown at the same times as the gas density images. 
Each figure displays the range in the phase as initial conditions (specifically, 
the gas and bulge fraction of the simulated M5 1 galaxy and orbital inclination 
of the satellite) are varied for a given concentration-scale radius value. The 
phase of the data of M5 1 is overplotted in each figure as a dot-dashed green 
line. 



We have investigated the dependence of the phase of the 
m = 1 mode at the best-fit time on parameters aside from 
the scale radius. Figure H] (a-c) depicts the variation (the 
spread is shaded in yellow) of the phase of the m = 1 mode 
for a given concentration-scale radius value as the bulge and 
gas fraction of the simulated M51 galaxy are varied, along 
with the orbital inclination of the satellite, for the three scale 
radius-concentration cases we have considered, namely (a) 
Rs=n,C=l4, to (b) Rs =n,C = 9.4, to (c) = 32,C = 5. 
Specifically, we consider here the fiducial model from CBCB 
that does not include a bulge, as well as a bulge mass frac- 
tion mt, = Mbuige/A^totai = 0.012. We consider gas fractions 
fg = 'Mgas/'Mdisk (where mdisk is the disk mass) of fg = 0.142 
(the fiducial case) and fg = 0.25. We consider the fiducial 
case from CBCB of a co-planar orbit for M5 1 's satellite, as 
well as an inclined orbit (the "e" orbit as denoted in CBCB). 

As is clear, while varying initial conditions has some effect 
on the phase, it does not change the essential character of the 
phase to the point that it would preclude us from discriminat- 
ing from the different scale radius cases. In each panel, the 
phase of the m = 1 mode of the HI data of M5 1 is overplotted 
as a dot-dashed green line for comparison. We found earlier 
in CBCB that varying the equation of state of the gas in the 
primary galaxy, and modestly varying the orbital velocity of 
the satellite, does not significantly affect our results, and we 
do not vary those parameters here. In the future, we will carry 
out a larger parameter survey and study the detailed effects of 
these and other parameters on this metric. Our initial study 
here does show that this metric is reasonably robust to the 
variation of other parameters, i.e., the shape of the gradient is 
sufficiently distinct to allow us to characterize the scale radius 
of the dark matter halo. 

We have shown the phase behavior here at the time when 
a given simulation best fits the Fourier amplitudes of the HI 
data. CB09 and CBCB have outlined how the best-fit simu- 
lation from a simulation parameter survey can be displayed 
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on a variance vs variance plot. Disturbances in the gas disk 
dissipate on the order of a dynamical time. As such, even 
the best-fit simulation traces a trajectory on the variance vs 
variance plot, with the minimum yielding the best-fit time. 
We have examined the time variation of these simulations as 
well. We do not find that the time variation of the phase would 
allow distinct scale radius cases to mimic each other. We con- 
clude therefore that the phase of the m= 1 mode allows us to 
discriminate between different scale radii within a factor of 
~ 1.5 (i.e., we can clearly discriminate the Rs = 17 case from 
the/?, = 11 case). 

3. DISCUSSION 

There are some caveats in our work here that are worth not- 
ing for the appUcation of this method in generality. We have 
not investigated in this paper the effect of non-spherical halos 
on either the Fourier amplitudes or the phase. Recent work 
has studied the effects of prolate halos (Banerjee & Jog 20 11) 
on the vertical structure of the Milky Way's HI disk. We defer 
a detailed investigation of these effects to a future paper. Pre- 
Uminary work (Chakrabarti, Debattista & Blitz, in prepara- 
tion) finds that the evolution of the halo shape (and thereby the 
strength of the Fourier amplitudes in the outskirts of the gas 
disk as determined by intrinsic processes) is significantly af- 
fected by gas cooling and angular momentum transport from 
the gas to the halo. Debattista et al. (2008) found earlier that 
the presence of a baryonic (stellar) disk affected the shape of 
the dark matter halo, leading to a rounder halo as the disk 
is grown. We have recently examined similar simulations 
which include a gaseous component (Chakrabarti, Debattista 
& Blitz, in prep) and find that the presence of the gaseous 
component renders halo shapes more spherical over time after 
the onset of gas cooling. The shapes of halos are considerably 
rounder close to present day, when the stellar to dark mat- 
ter masses are comparable to local spirals (M*/Mdm ~ 0.03) 
(Leroy et al. 2008), relative to j ~ 2. Examination of the 
Fourier amplitudes of these simulations close to present day 
indicates that the Fourier amplitudes in the outskirts are < 10 
%, which suggests that non-spherical halos would not be the 
primary contributor to the observed strength of the Fourier 
amplitudes in local spirals. This is certainly valid for M51, 
which is known to be a tidally interacting system (CBCB; 
Salo & Laurikainen 2000). Local spirals that display small (< 
10 % in the Fourier amplitudes, as defined in CB09, relative 
to the axisymmetric mode) perturbations in the outskirts, and 
have no visible tidally interacting companion may require a 
more careful treatment of the effect of halo shapes. Nonethe- 
less, the presence of a cold disssipative gaseous component 
needs to be modeled to more comprehensively understand the 
evolution of halo shapes. While cosmological simulations 
(Maccio et al. 2008) find halo shapes rendered non-spherical 
(in many cases triaxial) by repeated mergers, this effect may 
be mitigated in reality by a cold gaseous component trans- 
ferring its abundance of angular momentum to the dark halo, 
thereby rendering it more spherical. 

Our observationally motivated inference of the scale radius 
of the dark matter halo may be useful for a number of reasons. 
Firstly, the NFW density profile has some dependence on 
the spectral index of primordial density fluctuations (NFW96; 
NFW97). Secondly, this profile emerges within the context 
of dissipationless, cosmological simulations, and baryonic ef- 
fects may well serve to alter this profile. Therefore, an inde- 
pendent observationaUy motivated probe of the density profile 
for specific spiral galaxies can serve to support the theoreti- 



cal bases of this so-called universal density profile, wherein 
galaxies have grown hierarchically from primordial density 
fluctuations. Performing such a study for a large sample of 
spirals in the Local Volume, as can be done in the near future 
with galaxies from the THINGS survey, will allow for a statis- 
tical determination of the scale radii of spiral galaxies. Such 
a study can be performed to zeroth-order using the scaling re- 
lations in Chang & Chakrabarti (20 11) that allow one to infer 
the satellite mass from the observed HI map, and refined with 
numerical calculation. Another important tracer of the mass 
distribution is the observed rotation curve (Bosma 1978; Ru- 
bin et al. 1982; Blitz 1979). The rotation curve for many 
spirals cannot be determined beyond ~ 10 disk scale lengths. 
We leave a comparison of this method to the rotation curve to 
a future paper. It would also be worthwhile to compare our 
approach to gravitational lensing (Treu & Koopmans 2002; 
Mandelbaum et al. 2006), which also serves to probe mass 
distributions independent of the stellar light. Our method is 
ideally suited for analysis of disturbances in the outskirts of 
spiral galaxies in the Local Volume, and can in the future be 
appUed to a higher redshift population with observations from 
instruments such as SKA. 



4. CONCLUSION 

• We find that the shape of disturbances in the fragile, ex- 
tended HI disks of galaxies reflects the underlying density 
profile of the dark matter halo in spiral galaxies. We employ 
the phase of the m = 1 mode as our primary metric of com- 
parison to observed HI data. Building on prior results from 
CBCB where the satellite mass and pericentric approach dis- 
tance was determined for M51, we find that the scale radius 
of the dark matter halo can be determined from the phase of 
the m = 1 mode. 

• The three regimes of interest for the phase of the m = 1 
mode that delineate the variation of the scale radius and how 
it is reflected in the phase are: d/dr{(/)i) < for r/R^ < 1, a 
transitional region air/Rs^ 1, and d/dr((f)i) > for r/Rs > 1 
til the edge of the gas disk. Long baseUne HI observations 
that can reveal this transitional region between the shal- 
low profile (that would produce more loosely wrapped spirals) 
and the steeper profile (that would produce more tightly 
wrapped spirals) will allow us to quantitatively infer the scale 
radius of dark matter halos of specific spiral galaxies. The 
slight uncertainty in the phase that arises due to variations in 
initial conditions and orbits does not preclude us from infer- 
ring the value of the scale radius to within a factor of ^ 1.5. 
We have discussed preliminary work that suggests that other 
effects (such as intrinsic processes arising from non-spherical 
halos) will not significantly affect our results here. 

• Application of this method to the Whirpool Galaxy yields 
a scale radius of 17 kpc, which is consistent with the expected 
range of scale radii of ~ IO^^Mq halos in dissipationless cos- 
mological simulations. 

• We have also demonstrated that for r/R^ > 1, our results 
depend essentially on the scale radius and not on the concen- 
tration parameter. These results suggest that we can utilize 
observed disturbances in the extended HI disks of galaxies 
to characterize the density profile of the dark matter halo in 
spiral galaxies. In the future, we will apply this method to a 
large sample of spirals to obtain a statistical determination of 
the scale radii of spiral galaxies in the Local Volume. 
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